AD-A1U  791  VIRGINIA  POLYTECHNIC  INST  AND  STATE  UNIV  BLACKSBUR6  —ETC  F/6  20/A 

COMPUTATION  OF  HYPERSONIC  LAMINAR  VISCOUS  FLOW  OVER  A  BOOY  WITH— ETC<U> 
MAY  81  MO  KIM*  C  H  LEWIS 


UNCLASSIFIED 


NL 


ADA111791 


COMPUTATION  OF  HYPERSONIC  LAMINAR  VISCOUS  FLOW  OVER  A  BODY 
WITH  MASS  TRANSFER  AND/OR  SPIN  AT  ANGLE  OF  ATTACK 


*  ** 
M.  D.  KIM  and  C.  H.  LEWIS 


Aerospace  and  Ocean  Engineering  Department 


Virginia  Polytechnic  Institute  and  State  University 


Blacksburg,  Virginia  24061 

//  ft 

ABSTRACT 


Computational  results  of  hypersonic  laminar  viscous  flow  over 
blunt-nosed  bodies  with  mass  transfer  and/or  spin  are  presented.  An 
implicit-iterative  numerical  scheme  at  each  marching  step  is  used  to 
solve  the  parabolized  Navier-Stokes  equations.  The  code  takes  into 
account  the  mass-transfer  and  spin  effects  in  the  boundary  condition 
at  the  body  surface  and  the  periodic  effect  in  the  boundary  condition 
around  the  body.  To  facilitate  convergence,  the  v-momentum  equation 
has  been  used  rather  than  the  continuity  equation  to  provide  the 
pressure  at  the  body  boundary.  The  coupling  effects  of  mass  trans¬ 
fer  and  spin  on  the  integrated  forces  and  moments  as  well  as  the  sur¬ 
face  properties  are  presented  and  discussed. 


INTRODUCTION 


In  recent  years  the  analysis  of  hypersonic  viscous  flow  past  reentry  ve¬ 
hicles  with  mass  transfer  and/or  spin  has  been  a  problem  of  considerable 
interest  in  ballistics  because  of  its  relevance  in  slowly  spinning  planetary 
entry.  Previously  a  numerical  method  (VSL3D)1  based  on  a  viscous  shock-layer 
approach  was  applied  to  the  computation  of  heat-  and  mass-transfer  effects  on 
three-dimensional  viscous  shock-layer  flows  by  Murray  and  Lewis.2  The  VSL3D 
method  was  found  to  give  accurate  predictions  even  at  very  high  altitude  be¬ 
cause  the  code  takes  into  account  wall  and  shock  slip  effects.  Another  nu¬ 
merical  analysis  for  the  mass-transfer  effect  using  the  viscous  shock-layer 
theory  was  given  by  Whitehead  and  Davis.3  The  VSL3D  method  has  a  restriction 
on  the  angle  of  attack  for  the  entire  solution  around  the  body.  In  other 
words,  since  the  numerical  scheme  used  is  parabolic  in  the  circumferential 
direction,  it  cannot  treat  the  crossflow  separated  region.  In  the  above 
studies,  spin  effects  on  the  flowfield  were  not  treated. 


Recently  a  numerical  approach  for  the  hypersonic  viscous  flowfield  solu¬ 
tion  was  developed  by  Lubard  and  Helliwell1*  using  the  parabolized  Navier- 
Stokes  equations  (PNS).  The  PNS  method  has  been  applied  to  the  various 
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problems  predicting  hypersonic  viscous  flows  over  bodies  at  high  angles  of  at¬ 
tack.  For  instance,  Waskiewicz  and  Lewis5  applied  the  method  to  a  variety 
of  freestream  conditions  and  investigated  the  limit  and  applicability  of  the 
method.  They  found  that  the  PNS  method  could  produce  accurate  predictions  of 
the  flowfield  during  reentry  up  to  at  least  38-deg  angle  of  attack.  They  also 
discussed  the  following  effects  of  axial  pressure  gradient  models.  For 
problems  involving  high  Mach  number  and  Reynolds  numbers,  the  implicit  gradient 
model  provided  unique  solutions  regardless  of  the  stepsize  with  no  noticeable 
increase  in  the  number  of  iterations  required  for  convergence.  Problems  in¬ 
volving  large  adverse  streamwise  pressure  gradients  due  to  high  angle  of  at¬ 
tack,  large  mass  transfer,  or  small  Mach  and  Reynolds  numbers  can  cause  dif¬ 
ficulties  in  convergence  for  the  implicit  pressure  gradient  model.  In  such 
cases  an  explicit  or  zero  formulation  was  needed,  and  the  accuracy  of  the  so¬ 
lution  was  dependent  on  the  streamwise  stepsize.  In  their  work,  the  VPI  VSL3D 
code  was  satisfactory  for  the  preparation  of  an  initial  plane  data  which  was 
used  to  start  the  PNS  solution.  In  other  works,  Gogineni  and  Lewis6  reduced 
the  computing  times  required  for  an  entire  solution  by  combining  the  VSL3D 
and  PNS  methods  and  also  by  using  a  fast  implicit-iterative  technique  known  as 
the  Pseudo  Elimination  Method.  The  PNS  method  has  been  further  developed  by 
Helliwell  et  al.7*8  by  introducing  a  more  general  body-generator  coordinate 
system,  thus  extending  the  applicability  to  more  complex  bodies.  A  stream¬ 
line  coordinate  system  also  has  been  incorporated  in  the  PNS  method  by  Kim 
and  Lewis9  in  an  effort  to  extend  the  capability  to  treat  extremely  high 
angles  of  attack. 

Spin  effects  on  the  hypersonic  viscous  flowfields  past  conical  bodies 
have  been  computed  and  analyzed  by  Agarwal  and  Rakich.10  A  PNS  code  has  been 
modified  to  include  the  periodic  condition  for  the  flow  profiles  in  the  wind¬ 
ward  plane  to  account  for  the  asymmetry  in  the  flow  due  to  spin.  They  also 
introduced  a  new  criterion  for  defining  crossflow  separation  behind  spinning 
bodies  which  generalizes  the  Moore-Rott-Se?  s  criterion  for  two-dimensional 
unsteady  separation.  The  onset  of  separation  in  the  flowfield  is  characterized 
by  the  condition  3w/3y  =  0  at  w  =  0  where  w  is  the  circumferential  component 
of  velocity  and  y  the  coordinate  normal  to  the  body.  In  their  work,  detail 
analyses  were  made  concerning  the  effects  of  angle  of  attack,  freestream 
Reynolds  number,  nose  bluntness,  and  finally  spin  rate  on  the  Mngnus  force 
components  and  various  flow  properties.  Coning  motion  also  was  considered  by 
Agarwal  and  Rakich11  using  a  PNS  method  further  modified.  Another  study  of 
the  spin  effects  for  cones  or  ogive-cylinder-boattailed  bodies  has  been  re¬ 
ported  by  Sturek  and  Schiff12  using  another  PNS  approach.  None  of  the  fore¬ 
going  investigations  considered  the  coupling  effects  of  mass  transfer  and  spin. 

In  the  present  study,  computational  results  and  discussion  are  presented 
for  the  parametric  effects  of  mass  transfer  and/or  spin  together  with  the 
reasonable  computing  times  consumed  for  each  test  case.  The  original  PNS 
code  used  the  continuity  equation  to  obtain  the  wall  pressure  at  the  body 
boundary  in  the  case  of  mass  transfer.  Theoretically  this  method  appears 
quite  reasonable  because  it  can  provide  a  finite  value  of  the  wall  pressure 
when  the  normal  velocity  at  wall  <}°es  not  vanish,  and  moreover  it  is  a 
relatively  simple  equation.  However,  by  numerical  experiment  it  has  been 
found  that  a  solution  is  hardly  obtainable  by  this  method  in  the  case  of  mass 
transfer.  Therefore,  in  the  present  work,  the  v-momentum  equation  is  dif¬ 
ferenced  to  provide  the  wall  pressure  at  the  body  boundary.  In  most  cases. 
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the  flowfield  solution  has  been  obtained  without  convergence  problems  by  the 
current  method. 

In  the  next  section,  the  theoretical  background  and  procedure  will  be  pre¬ 
sented  which  includes  a  discussion  of  the  governing  equations,  coordinate  sys¬ 
tem,  boundary  conditions,  and  numerical  solution.  A  few  methods  to  provide 
initial  plane  data  are  also  introduced,  and  the  effect  of  initial  plane  data 
for  a  starting  solution  is  discussed.  In  the  present  work,  the  zero  pressure 
gradient  model  has  been  used  for  all  the  test  cases  in  order  to  obtain  a  so¬ 
lution  in  the  case  of  high  angles  of  attack  with  mass  transfer  and/or  spin. 

Two  sets  of  freestream  conditions  are  considered.  The  first  case  is 
Mach  =  5,  a  =  2  and  Re„j/ft  =  6.21  x  10^.  The  second  case  is  Mach  =  18,  a  =  15 
and  Reoo/ft  =  5.45  x  10^.  For  both  cases,  the  body  geometry  is  a  spherically 
blunted  7-deg  half-angle  cone  0.817-ft  long.  Spin  rates  considered  are  2000 
and  8000  rev/min  and  mass-transfer  rates  0.00125  and  0.00250.  For  Case  1,  an 
asymmetric  sinusoidal  distribution  of  mass  transfer  around  body  as  well  as  a 
constant  distribution  has  been  considered.  For  each  freestream  condition, 
computational  results  were  obtained  for  four  parametric  cases,  i.e.,  basic 
case,  spin  only  case,  mass  transfer  only  case,  and  finally  mass  transfer  and 
spin  case.  Parametric  comparisons  are  presented  to  study  the  effect  of  mass 
transfer  and/or  spin  on  the  Magnus  force  components  as  well  as  all  the  surface 
properties.  It  will  be  noted  that  a  significant  coupling  effect  on  the  Magnus 
force  can  result  from  mass  transfer  and  spin.  Reasonable  computing  times  re¬ 
quired  for  the  solutions  are  also  presented. 

ANALYSIS 

GOVERNING  EQUATIONS 

The  governing  equations  are  the  laminar,  steady  parabolized  Navier-Stokes 
equations  derived  by  Lubard  and  Helliwell4  under  the  assumption  that  the 
gradients  of  the  shear  stress  in  the  streamwise  direction  are  negligibly  small 
compared  with  the  gradients  in  the  direction  of  normal  and  circumferential  di¬ 
rections.  The  resulting  equations  are  parabolic  in  the  streamwise  direction 
and  elliptic  in  the  crossflow  direction.  The  governing  equations  include  the 
continuity,  three  momentum,  and  energy  equations.  The  fluid  density  is  given 
by  the  equation  of  state  for  a  perfect  gas.  Prandtl  number  and  specific  heat 
are  assumed  constant  and  Sutherland's  viscosity  law  is  usee 

In  addition  to  the  usual  boundary-layer  order  of  magn.'r. ....  terms,  d<J>  is 
assumed  to  be  0(e)  to  include  the  crossflow  separation  and  the  i ormal  velocity 
v  is  taken  as  0(1).  The  flow  variables  are  nondimensionalized  by  their 
freestream  values.  The  normal  coordinate  is  normalized  by  the  shock-layer 
thickness  to  facilitate  the  computation  of  the  shock  location.  The  complete 
set  of  governing  equations  together  with  the  derivation  procedure  can  be  found 
in  Reference  4. 

In  the  development  of  these  equations,  a  body-oriented  coordinate  system 
is  used  (Figure  1),  with  x  taken  along  the  body  generator,  y  normal  to  the 
body  surface,  and  $  the  circumferential  coordinate  where  =  0  represents 
the  windward  ray.  The  coordinate  system  does  not  rotate  with  the  body. 


BOUNDARY  CONDITIONS 


The  following  equations  are  used  for  the  boundary  conditions  at  the  body 
surface  with  mass  transfer  and  spin. 

u  =  0 

v  =  specified  distribution  (mass  transfer) 
w  =  specified  distribution  (spin) 
h  =  specified  distribution 
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To  obtain  the  wall  pressure  for  the  mass-transfer  and  spin  case,  the 
original  code  utilized  the  continuity  equation;  however,  the  method  turned 
out  to  be  inefficient  for  the  convergence  of  a  solution  by  a  numerical  ex¬ 
periment.  It  has  been  found  that  the  employment  of  the  v-momentum  equation 
gives  a  more  stable  solution.  Thus,  for  the  present  analysis  the  v-momentum 
equation  has  been  taken  and  differenced  using  a  one-sided  differencing  scheme 
to  provide  the  pressure  at  the  body  boundary. 

In  <  rder  to  obtain  the-  outer  boundary  conditions,  Rankine-Hugoniot  jump 
conditions  are  utilized  at  the  shock.  The  components  of  freestream  velocity 
in  the  body  oriented  coordinate  system  are: 

u  =  cos  0  cos  a  -  sin  0  sin  a  cos  4> 
c  c 

v  =  -sin  0  cos  a  -  cos  6  sin  a  cos  6 
c  c 

w  =  sin  <J>  sin  a 

From  the  freestream  components  and  jump  conditions,  five  conservation  equa¬ 
tions  are  obtained  which  can  be  used  to  determine  the  aftershock  properties. 

To  uniquely  determine  the  six  unknowns  £,  u,  v,  w,  p,  h  at  the  shock  boundary, 
one  more  equation  is  required.  Thus,  the  one-sided  differencing  of  the 
continuity  equation  provides  the  sixth  equation. 

Since  the  windward  and  leeward  surfaces  are  not  symmetry  planes  for  the 
spinning  cone,  a  periodic  condition  for  the  flow  profiles  in  the  windward 
plane  is  specified,  thus  providing  the  boundary  conditions  in  the  ^-coordinate 
direction. 
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In  a  rectangular  y-4>  grid,  the  shock  may  not  fall  on  a  mesh  point.  To 
resolve  this  difficulty,  a  coordinate  transformation,  ri  =  y /£  (x,4>)  is  made, 
where  £(x,<j>)  is  local  shock-layer  thickness.  Due  to  the  coordinate  trans¬ 
formation,  the  shock-layer  thickness  E,  appears  in  all  the  governing  equa¬ 
tions,  thus  requiring  one  more  governing  equation,  3£/3 y  =  0  to  be  differenced. 

INITIAL  CONDITIONS 

For  a  numerical  flowfield  solution  which  utilizes  a  marching  scheme, 
preparation  of  an  accurate  initial  data  plane  (IDP)  is  one  of  the  most 
crucial  conditions  for  a  successful  start  of  solution.  By  previous  investi¬ 
gations5’6  the  viscous  shock-layer  method  (VSL3D)1  for  blunt  bodies  was  found 
to  be  able  to  generate  a  satisfactory  initial  data  plane  to  start  the  PNS  so¬ 
lution.  Thus,  the  entire  flow  properties  including  the  shock  shape  must  be 
supplied  at  an  initial  data  plane  to  get  the  PNS  code  started. 

Three  different  methods  have  been  used  to  generate  an  initial  data  plane. 
The  first  method  is  to  construct  a  one-step  IDP  using  the  VSL3D  method  at  the 
sphere-cone  juncture.  In  most  cases,  the  one-step  IDP  constructed  by  the 
VSL3D  code  was  satisfactory  to  start  the  PNS  solution.  The  one-step  IDP  is 
obtained  by  an  axisymmetric  VSL  solution  and  necessary  rotations  of  the  solu¬ 
tion.  In  this  case,  a  few  initial  marching  stepsizes  should  be  small  to  start 
a  solution. 


The  second  method  is  to  prepare  a  two-step  IDP  using  VSL3D  code  at  any 
desired  axial  station  of  the  body.  This  method  is  better  than  the  first  one 
for  the  convergence  of  the  starting  solutions,  and  hence  the  initial  marching 
step-sizes  can  be  taken  larger.  The  third  method  is  to  prepare  a  two-step 
IDP  using  LUB2  code  which  is  another  PNS  code  developed  by  Agopian  et  al.7  for 
more  general  bodies.  In  this  method,  an  axisymmetric  VSL  solution  provides 
an  IDP  at  the  blunt  nose  for  the  LUB2  solution,  and  then  the  LUB2  code  again 
prepares  an  IDP  at  a  desired  axial  station  by  the  necessary  three-dimensional 
calculations.  This  method  is  somewhat  tedious  but  can  produce  an  accurate  IDP 
even  when  the  angle  of  attack  is  over  30  degrees.  In  the  preparation  of  IDP 
for  the  test  cases,  the  spherical  nose  tips  of  the  bodies  were  considered  to 
have  neither  spin  nor  mass  transfer. 

NUMERICAL  SOLUTION 


The  equations  are  solved  by  implicit  differencing  in  the  y,  (£-plane.  The 
x  derivatives  are  approximated  by  a  backward  difference  while  y  and  <}> 
derivatives  use  an  unequally  spaced  three-point  difference  formula,  e.g.. 


3f 

dx 


and 


v. 


3f 

3y  *  Vl  fk-l  +  akfk.  +  ak+l  fk+l 


11-195 


- '  ion  Fop 

;  i  '  CR'i&I 
t  '  3 

'pd 

1  •’  cnt  1  oti_ 

n 

□ 

.  ■ ut i on/ 

-  i 

. . .  _.j 

•  •  ■  v-ibii  5 \y 

Codes 

A  -’ft  ‘  1  flu 

>  /  *- 

;  Sm'r" 

and 


where  the  are  determined  from  the  Taylor  series  expansions  of  f^ 
f^+^  around  the  k-th  location,  and  hence  they  are  functions  of  the  stepsize 
between  the  grid  points. 

The  PNS  equations  and  the  perturbation  equations  during  the  iterations 
can  be  written  in  the  matrix  form 

F(U)  =  0 

An  =  Un+1  -  Un 

where  F  denotes  the  governing  equations  for  the  unknowns  (J,  and  is  the 

column  matrix  of  perturbation  properties  at  n-th  iteration. 

After  differencing,  the  equations  are  linearized  by  the  Newton-Raphson 
method.  Thus,  all  the  nonlinear  terms  of  perturbation  properties  are  dropped. 
This  results  in  an  equation  of  the  form 

M  (Un)  An  =  -  F(Un)  ;  n  =  0,  1,  2,  ... 

with  an  initial  guess  to  the  solution  of  the  governing  equations,  and 
M  the  Jacobian  of  f\  Because  of  the  size  of  the  system,  instead  of  solving 
the  above  equation  directly,  the  Gauss-Seidel  iteration  method  is  used.  In 
the  ^-direction,  an  implicit-iterative  scheme8  is  used,  but  after  convergence, 
the  solution  obtained  is  a  fully  implicit  one.  After  convergence,  another 
step  in  the  marching  direction  is  taken  and  the  whole  procedure  is  repeated. 

Convergence  problems  often  occur  when  the  marching  solution  encounters 
a  strongly  separated  region  due  to  crossflow.  To  reduce  the  problem,  either 
zero  or  explicit  pressure  gradient  model  may  be  used.  Smaller  marching  step- 
sizes  and/or  pressure  smoothing  in  the  y  coordinate  direction  also  can  help 
the  convergence.  When  the  angle  of  attack  is  very  high,  more  careful  selection 
of  grid  sizes  and  preparation  of  more  accurate  initial  plane  data  are  required 
in  order  to  start  a  solution. 

Mass  transfer  and  spin  add  more  disturbances  in  the  viscous  flowfield, 
thus  deteriorating  the  convergence  of  the  solution.  A  large  mass-transfer 
distribution  along  the  body  can  cause  an  axial  flow  separation  at  some  down¬ 
stream  station  which  causes  the  marching  solution  to  fail.  To  reduce  the 
marching  stepsize  problem  in  a  solution,  the  marching  stepsizes  are  controlled 
internally  by  the  code  considering  the  number  of  iterations  taken  for  the 
solution  at  the  previous  step. 

RESULTS  AND  DISCUSSION 

For  the  present  study,  a  body-generator  coordinate  system  is  used  with 
x  along  the  body,  y  normal  to  the  body  and  <f>  around  the  body,  and  it  does  not 
rotate  with  the  body.  Figure  1  shows  the  coordinate  system  and  the  sign 
conventions  for  aerodynamic  coefficients.  Since  this  sign  convention  is  not  a 
universal  one,  it  is  recommended  to  be  used  only  for  the  interpretation  of 
the  directions  of  the  forces  and  moments  reported  in  this  paper. 


The  initial  plane  data  for  the  test  cases  have  been  obtained  using  a  per¬ 
fect  gas  viscous  shock-layer  approach  (VSL)  at  the  sphere-cone  tangency  point, 
and  all  the  flow  properties  together  with  forces  and  moments  are  transferred 
to  the  PNS  solution.  To  start  the  PNS  solution,  an  initial  marching  step- 
size  of  less  than  one-tenth  of  the  nose  radius  has  been  used,  and  it  was  in¬ 
creased  internally  in  the  code,  considering  the  number  of  iterations  taken 
and  the  maximum  permissible  marching  stepsize  for  a  meaningful  solution.  In 
the  present  calculations,  101  grid-points  in  the  normal  direction  and  19 
planes  around  the  body  are  used  for  all  cases. 

The  original  PNS  code  utilized  the  continuity  equation  in  order  to  pre¬ 
dict  the  wall  pressure  at  the  body  boundary  when  the  body  had  mass  transfer. 
However,  by  numerical  experiment,  it  has  been  found  that  the  method  is  likely 
to  give  an  unstable  and  divergent  solution  even  in  the  case  of  very  small 
mass  transfer.  Therefore,  in  the  present  work  the  v-momentum  equation  is 
taken  and  differenced  at  the  body  boundary.  By  this  method,  stable  solutions 
have  been  obtained  in  most  cases  with  moderate  mass  transfer.  A  large  mass 
transfer  at  the  wall  can  cause  an  axial  flow  separation  at  some  downstream 
station.  In  such  a  case,  it  should  be  noted  that  the  PNS  method  cannot  treat 
the  primary  flow  separated  region. 

In  the  case  of  high  angle  of  attack  with  large  mass  transfer  and  spin, 
a  relatively  large  amount  of  perturbations  can  occur  in  the  viscous  flowfield. 
These  perturbations  have  an  adverse  effect  on  the  convergence  of  a  solution. 

For  this  reason,  in  the  present  computations,  the  zero  pressure  gradient  model 
in  the  marching  direction  was  used  to  improve  the  convergence.  The  effect  of 
the  three  different  kinds  of  axial  pressure  gradient  models  has  been  dis¬ 
cussed  in  detail  by  Waskiewicz  and  Lewis.7 

To  investigate  the  effect  of  mass  transfer  and/or  spin  on  the  viscous 
flow  past  a  sphere-cone  geometry,  two  sets  of  freestream  conditions  have  been 
chosen  and  parametric  computations  have  been  made.  Case  1  considers  a  free¬ 
stream  condition  of  M  =  5,  a  =  2  and  Re^/ft  =  6.21  x  105.  Case  2  considers 
a  freestream  condition  of  M  =  18,  a  =  15  and  ReaJtt  =  5.45  x  10s.  The  body 
geometry  is  a  spherically  blunted  7-deg  half-angle  cone  0.817  ft  long.  For 
both  cases,  a  parametric  comparison  is  presented  to  study  the  effect  of  spin, 
the  effect  of  mass  transfer,  and  the  coupling  effect  of  mass  transfer  and 
spin.  For  Case  1,  the  effect  of  the  asymmetric  mass-transfer  distribution 
around  the  body  as  well  as  constant  mass-transfer  distribution  will  be 
analyzed  in  detail.  Spin  rates  considered  are  2000  and  8000  rev/min.  Case 
la  has  a  sinusoidal  mass-transfer  distribution  around  the  body,  and  the 
maximum  value  of  mass-transfer  rate  is  0.0025.  Case  lb  and  Case  2  have  a 
mass  transfer  rate  of  0.00125,  which  is  constant  along  and  around  the  body. 
Details  of  the  freestream  conditions  are  given  in  Table  1. 

For  Case  1,  the  flowfield  solution  over  the  body  with  mass  transfer 
could  not  be  obtained  up  to  the  body  end.  For  Case  la  with  mass  transfer  only, 
the  solution  was  obtained  up  to  s  =  24.8  where  axial  flow  separation  occurred 
due  to  the  mass-transfer  effect.  For  Case  la  with  both  mass  transfer  and  spin, 
the  solution  was  obtained  up  to  s  *  23.3  before  axial  flow  separation  occurred. 
The  effect  of  mass  transfer  on  axial  flow  separation  will  be  further  dis¬ 
cussed  later. 


Figure  2  shows  a  sinusoidal  distribution  of  mass  transfer  for  Case  la, 
which  is  designed  to  exert  a  relatively  large  side  force  in  the  minus  direc¬ 
tion.  In  Figure  3,  the  Magnus  force  due  to  the  sinusoidal  mass-transfer 
distribution  is  shown,  and  also  it  can  be  observed  that  the  coupling  effect 
of  mass  transfer  and  spin  on  the  Magnus  force  is  remarkable.  Figure  4  shows 
the  Magnus  force  components  distribution  along  the  body  for  Case  la,  where  it 
is  noted  that  the  Magnus  force  component  due  to  the  wall  pressure  distribution 
constitutes  the  largest  part  of  the  total  Magnus  force,  while  the  components 
due  to  the  axial  and  the  circumferential  shear  stresses  are  negligible. 

For  Case  lb,  we  have  constant  mass  transfer  around  the  body.  Hence  in 
this  case,  the  Magnus  force  due  to  mass  transfer  should  be  negligible.  In 
Figure  5,  it  is  noted  that  the  Magnus  force  due  to  the  coupling  effect  of 
the  mass  transfer  and  spin  is  substantial,  and  the  effect  is  over  three  times 
as  much  as  the  spin-only  effect  on  the  Magnus  force.  The  Magnus  force  com¬ 
ponent  due  to  the  wall  pressure  prevails  over  the  others  due  to  the  wall  shear 
stresses  for  this  case  as  shown  in  Figure  6. 

Figure  7  shows  the  effect  of  mass  transfer  and  spin  on  the  Magnus  force 
for  Case  2  which  has  relatively  large  Mach  number  and  high  angle  of  attack. 

For  this  case,  it  is  noted  that  the  coupling  effect  of  mass  transfer  and  spin 
is  again  remarkable.  Therefore,  on  the  basis  of  the  present  computational 
results,  it  is  concluded  that  this  coupling  effect  should  not  be  neglected 
in  a  flight  dynamics  analysis  of  slow-spin  reentry  vehicle  with  mass  transfer. 
Figure  8  shows  that  the  Magnus  force  component  due  to  the  crossflow  shear 
stress  is  not  negligible,  and  it  has  a  different  direction  from  the  F_  for 
the  Case  2.  "w 

The  effect  of  mass  transfer  and  spin  on  the  wall  pressure  is  shown  in 
Figure  9,  where  the  spin  effect  was  almost  negligible.  Hence,  the  differences 
shown  are  largely  due  to  mass  transfer.  It  is  noted  that  the  mass-transfer 
effect  increases  the  wall  pressure  because  the  thickening  of  the  viscous  layer 
increases  the  displacement-induced  pressure  interaction.  In  this  figure,  an 
unstable  solution  due  to  axial  flow  separation  can  be  seen  at  s  =  24  and 
<}>  =  180  deg. 

Figure  10  shows  the  effect  of  asymmetric  mass  transfer  and  spin  on  the 
wall  pressure  at  s  =  23.  In  fact,  the  pressure  difference  between  the  two 
cases  is  due  mostly  to  the  mass-transfer  effect.  Since  Case  la  has  maximum 
mass  transfer  in  the  <p  =  270-deg  plane  of  the  body,  the  pressure  difference 
is  larger  on  this  plane  of  the  body,  as  can  be  seen  in  the  figure.  Figure  11 
shows  the  effect  of  mas^  transfer  and  spin  on  the  shock-layer  thickness  which 
slightly  increases  the  shock-layer  thickness.  Figure  12  shows  the  effect  of 
mass  transfer  and  spin  on  the  shock-layer  thickness  around  the  body  at  s  =  23 
where  the  slight  increase  in  the  shock-layer  thickness  is  due  mostly  to  the 
mass  transfer  effect,  not  the  spin  effect. 

The  effect  of  mass  transfer  and  spin  on  the  surface  heat  transfer  along 
the  body  is  shown  in  Figure  13.  Mass-transfer  effect  decreases  the  wall 
heat  transfer,  and  the  effect  is  more  sizable  on  the  windward  side.  Figure  14 
shows  the  effect  of  asymmetric  mass  transfer  and  spin  on  the  heat  transfer  at 
s  *  23  for  Case  la.  The  big  asymmetric  difference  between  the  two  cases  is 
largely  due  to  the  mass-transfer  effect.  Figures  15  and  16  show  the  effects 
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of  mass  transfer  and  spin  on  the  streamwise  wall  shear,  and  the  effects  are 
quite  similar  to  the  effects  on  the  heat  transfer.  Figure  17  shows  the  ef¬ 
fects  of  asymmetric  mass  transfer  and  spin  on  the  crossflow  wall  shear  at 
s  =  23.  The  effect  is  to  decrease  the  shear  stress  in  absolute  value  on  both 
sides  of  the  body.  As  previously  mentioned,  for  Case  1  the  spin  effects  on 
wall  pressure,  shock-layer  thickness,  heat  transfer,  and  wall  shear  were  very 
small  compared  to  the  mass-transfer  effects.  However,  for  Case  2  which  has  a 
larger  spin  rate  than  Case  1,  the  spin  effect  was  not  negligible,  which  will 
be  shown  later. 

Figure  18  shows  the  effect  of  mass  transfer  on  axial  flow  separation 
for  Case  la.  In  the  figure  it  can  be  observed  that  the  mass  transfer  retards 
the  axial  flow  velocity  in  the  boundary  layer.  The  effect  becomes  larger 
further  downstream,  and  finally  the  axial  flow  separates  at  s  =  24  as  can  be 
seen  in  the  figure.  The  PNS  method  cannot  treat  the  axial-flow  separated 
region,  hence  the  solution  diverged  at  this  station.  In  the  downstream  re¬ 
gion  before  axial  flow  separation,  the  numerical  solution  took  a  relatively 
large  number  of  iterations  and  small  marching  stepsizes.  This  effect  resulted 
in  the  relatively  large  computing  time  for  this  case  as  can  be  seen  in  Table  3. 

Case  2  had  a  spin  rate  of  8000  rpm  and  constant  mass-transfer  rate  of 
0.00125.  The  computational  results  for  this  case  showed  more  sizable  effects 
of  spin  on  the  wall  pressure  and  other  surface  properties.  Hence,  parametric 
comparisons  will  be  presented  for  this  case.  Figures  19  and  20  show  the  wall 
pressure  in  case  of  mass  transfer  and  spin.  Acutally  it  was  hard  to  distinguish 
each  effect  separately  in  this  kind  of  figure.  Hence,  it  was  desirable  to 
compare  the  ratios  of  each  pressure  to  its  basic  value.  Here,  the  basic  value 
means  the  computed  wall  pressure  for  the  body  with  neither  mass  transfer  nor 
spin.  Thus,  Figure  21  shows  the  coupling  effect  as  well  as  the  separate  ef¬ 
fects  of  mass  transfer  and/or  spin  on  the  wall  pressure.  In  this  figure,  the 
spin  effect  appears  antisymmetric,  while  the  mass-transfer  effect  is  symmetric 
around  the  body.  The  coupling  effect  is  quite  clear,  and  it  has  nearly  ad¬ 
ditive  characteristic.  The  spin  effect  decreases  the  surface  pressure  in 
<t>  =  0  -  180-deg  of  the  body,  and  increases  the  pressure  in  <j)  =  180  -  360-deg 
of  the  body,  thus  producing  a  negative  Magnus  force  as  can  be  expected.  The 
symmetric  mass  transfer  itself  cannot  produce  a  Magnus  force,  but  it  in¬ 
creases  the  wall  pressure  symmetrically.  From  Figures  20  and  21  it  is  ob¬ 
served  that  at  station  s  =  30.06  a  strongly  separated  flow  is  developed  in 
the  crossflow  direction.  For  this  case  at  ot  =  15-deg,  the  crossflow  separa¬ 
tion  begins  as  early  as  s  =  6,  as  can  be  seen  in  Figures  32  and  33. 

Figures  22  and  23  show  the  shock-layer  thicknesses  along  and  around  the 
body  respectively  where  the  effects  of  mass  transfer  and  spin  on  the  shock- 
layer  thickness  appear  negligible.  Figure  24  shows  the  effect  of  mass  trans¬ 
fer  and  spin  on  the  heat  transfer  in  the  windward  and  leeward  planes.  A  para¬ 
metric  comparison  of  the  effect  of  mass  transfer  and/or  spin  on  the  heat 
transfer  around  the  body  at  s  =  30.06  is  presented  in  Figure  26,  and  it  can 
be  interpreted  on  the  basis  of  the  data  given  in  Figure  25.  The  spin  effect 
on  the  heat  transfer  around  the  body  turns  out  quite  complex  due  to  the  cross- 
flow  separation  as  shown  in  Figure  26.  The  mass-transfer  effect  is  sym¬ 
metric  and  decreases  the  wall  heat  transfer.  The  coupling  effect  is  approxi¬ 
mately  additive  of  both  effects.  Figures  27  -  29  show  the  effects  of  mass 
transfer  and/or  spin  on  the  streamwise  wall  shear,  and  the  effects  have 
similar  trends  to  the  effects  on  the  heat  transfer. 
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Figures  30  and  31  show  the  effects  of  mass  transfer  and/or  spin  on  the 
crossflow  wall  shear  around  the  body  at  s  =  30.06.  The  effects  frequently 
change  their  directions  due  to  crossflow  separation.  The  coupling  effect 
is  again  additive  of  both  effects.  Since  the  current  PNS  method  is  elliptic 
in  the  crossflow  direction,  it  can  solve  the  crossflow  separated  region.  But 
the  crossflow  separation  produces  more  perturbations  in  the  flowfield,  and 
the  perturbations  have  an  adverse  effect  on  the  convergence  of  a  solution. 
Therefore,  generally  smaller  marching  stepsizes  are  required  to  pass  through 
the  region  where  the  crossflow  separation  begins.  Mass  transfer  causes  cross- 
flow  separation  at  an  earlier  station  on  the  body  while  spin  suppresses  the 
crossflow  separation.  A  definition  of  crossflow  separation  in  the  case  of 
spinning  body  has  been  introduced  by  Agarwal  and  Rakich.10  The  criterion 
generalizes  the  Moore-Rott-Sears  criterion  for  two-dimensional  unsteady  sepa¬ 
ration.  The  onset  of  separation  in  the  flowfield  is  characterized  by  the 
condition  3w/3y  =  0  at  w  =  0  where  w  is  the  circumferential  component  of 
velocity  and  y  is  the  coordinate  normal  to  the  body. 

Figure  32  shows  the  effect  of  spin  on  crossflow  separation  at  s  =  6.4 
and  <J>  =  160-deg  for  Case  2.  The  basic  case  already  has  the  separated  cross- 
flow  profile  at  this  station.  But  for  the  spin  case,  the  profile  is  not 
separated  yet  at  this  station  due  to  the  definition  introduced  above.  Figure 
33  shows  the  effect  of  mass  transfer  on  crossflow  separation.  It  is  observed 
that  mass  transfer  causes  crossflow  separation  at  an  earlier  station  on  the 
body.  As  previously  mentioned,  mass  transfer  can  cause  axial  flow  separation 
also.  In  a  numerical  experiment,  mass  transfer  of  a  few  percent  of  p  U  was 
found  to  easily  cause  an  axial  flow  separation  within  a  few  nose  radii  down¬ 
stream.  In  the  case  of  low  mass  transfer  as  in  Case  la,  the  axial  flow  was 
separated  at  s  =  24. 

Table  2  summarizes  all  the  force  and  moment  coefficients  obtained  from 
the  parametric  computations  for  the  two  test  cases.  The  effect  of  the  spin 
on  the  axial  and  normal  forces  is  found  to  be  negligible  while  the  effect  on 
the  Magnus  side  force  is  sizable.  An  asymmetric  mass  transfer  also  can  pro¬ 
duce  a  finite  amount  of  side  force.  From  the  table,  it  is  also  observed  that 
the  mass  transfer  decreases  the  normal  and  axial  forces.  An  important  fact 
found  in  the  present  work  is  that  mass  transfer  and  spin  have  a  measurable 
coupling  effect  on  the  Magnus  side  force,  which  is  shown  in  Table  2  as  well 
as  in  the  figures.  However,  no  experimental  data  are  known  to  the  authors 
at  this  time  to  validate  the  present  results  for  the  coupling  effects. 

Computing  times  for  these  cases  are  given  in  Table  3.  These  times  were 
obtained  on  the  IBM  3032  system  at  Virginia  Polytechnic  Institute  and  State 
University.  Most  cases  consumed  less  than  40  minutes  of  CPU  time  for  the 
entire  solution,  but  the  solutions  for  Case  la  and  lb  with  mass  transfer  took 
relatively  longer  computing  times  because  of  the  small  stepsizes  taken  before 
axial  flow  separation. 

SUMMARY 

The  viscous  flowfield  over  a  sphere-cone  undergoing  mass  transfer  and/or 
spin  at  angle  of  attack  has  been  calculated  using  the  parabolized  Navier- 
Stokes  equations.  In  the  mass-transfer  case,  the  wall  pressure  is  obtained 
using  the  v-momentum  equation  rather  than  continuity  equation  to  facilitate 


the  convergence.  The  main  results  of  the  present  work  can  be  summarized  as 
follows : 

1.  The  coupling  effect  of  mass  transfer  and  spin  on  the  Magnus  force 
is  remarkable,  i.e.,  it  is  much  more  than  the  additive  effect. 

2.  Mass  transfer  has  a  tendency  to  reduce  the  axial  and  normal  forces, 
and  it  increases  the  wall  pressure  and  decreases  the  heat  transfer  and  wall 
shear  stresses. 

3.  Mass  transfer  causes  earlier  axial  and  crossflow  separations. 

4.  Spinning  motion  has  negligibly  small  effect  on  the  normal  and 
axial  forces,  and  it  produces  asymmetric  effects  on  the  wall  pressure  as  well 
as  the  other  surface  properties  around  the  body. 

5.  Spin  has  a  tendency  to  suppress  the  appearance  of  crossflow  separa¬ 
tion  at  high  angles  of  attack. 
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NOMENCLATURE 


C 

P 

CFS 

CFW 


FY 

H 

h 

L 

M 


constant  pressure  specific  heat 

streamwise  skin-friction  coefficient,  2t*/(P  U?) 

'  ' r  oo  oo' 

*  2 

crossflow  skin-friction  coefficient,  2t^/ (p^U^) 
axial  force  coefficient 
normal  force  coefficient 
Magnus  force  coefficient 

pitching  moment  coefficient  about  blunt  nose-tip 

Magnus  moment  coefficient  about  blunt  nose-tip 
2  x  2  w 

sin  0  /  /  x  sin  <p  T  (x , <f>)  dd>  dx 

C  0  0 

4  sin20  P  P  x  sin  <J>  p  (x,4>)  d$  dx 
1  c  0  0  W 

sin  0  /X  P  x  cos  t,(x,4>)  d<f>  dx 

c  0  o  41 

total  Magnus  side  force,  F  +  F  +  F 

x  pw  T<f> 

2 

total  enthalpy,  H*/U 

* 

static  enthalpy,  h*/(C  T  ) 

p  00 

reference  length,  equivalent  cone  slantwise  length 
freestream  Mach  number 


11-202 
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-  1 


mass-transfer  rate  at  wall,  p*  v*/p  U 

w  w  00  00 

p  nondimensional  pressure,  P*/(P£JOU^) 

2 

Pw  wall  pressure,  pVp^ 

Pr  Prandtl  number 

3 

q  heat-transfer  rate,  <?*/(P00^00) 

r  local  body  radius,  r*/L 

R  dimensional  body  nose  radius  of  curvature 

n 

Re  /ft  freestream  unit  Reynolds  number 

oo 

ST1NF  Stanton  number,  q*/[p^U  (H*  -  H*)l 

s  surface  distance  coordinate  measured  along  the  body  from  nose 

tip,  s*/R 

n 

T  temperature,  T */Tga 

freestream  stagnation  temperature 

2  2  * 

T  ,  reference  temperature,  (  y~  1)  M  T  or  U  /C 
ref  r  oo  oo  p 

*  , 

T  wall  temperature,  T  / T 

w  w  oo 

U  dimensional  freestream  velocity 

OO 

u  velocity  in  the  x  direction,  u*/U 

v  velocity  in  the  y  direction,  v*/^ 

w  velocity  in  the  $  direction,  wVU^ 

x  coordinate  along  the  body  surface  with  constant  <j>,  x*/L 

y  coordinate  normal  to  the  surface,  y*/L 

YSH  bow-shock  standoff  distance  divided  by  R 

Z  center  of  pressure  from  blunt  nosetip  in  percentage  of  the  blunt 

rP  body  axial  length 

a  angle  of  attack,  degree 

Y  ratio  of  specific  heats 

q  transformed  normal  coordinate,  y/£ 


£ 

0 

C 

U 

K 

p 

T 


2  * 

perturbation  parameter,  £  =  y  (T^^)/ (P^U^lO 


cone  half-angle,  degree 


coefficient  of  viscosity,  y  /y 


bow-shock  standoff  distance  divided  by  L 


density,  p /p 


*  2 

primary  flow  wall  shear  stress,  t  / (p  U  ) 

X  oo  oo 


4> 


*  2 

crossflow  wall  shear  stress,  X  ,/ (p  U  ) 


<J>  circumferential  coordinate 

ft  angular  velocity  (rev/min) 

Subscript 

w  wall  value 

0  stagnation  condition 

°°  freestream  condition  (dimensional  quantity) 

Superscript 

*  dimensional  quantity 
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Computing  Times*  for  Test  Cases 
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Figure  5:  Case  lb  Effects  of  constant  mass  transfer  and  spin  on  the 
Magnus  force  along  the  body. 
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Figure  7:  Case  2  Effects  of  constant  mass  transfer  and  spin  on  the 
Magnus  force  along  the  body. 
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Figure  10:  Case  la  Effects  of  asymmetric  mass  transfer  and  spin  on  the 
surface  pressure  around  the  body  at  s  =  23. 
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Case  la  Effects  of  mass  transfer  and  spin  on  the  streamwise 
shear  along  the  body. 


T 


T 


T 


Figure  16:  Case  la  Effects  of  asymmetric  mass  transfer  and  spin  on 
the  streamwise  shear  around  the  body  at  s  *  23. 
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Figure  17:  Case  la  Effects  of  asymmetric  mass  transfer  and  spin  on 
the  crossflow  shear  around  the  body  at  s  =  23. 
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Figure  19:  Case  2  Effects  of  mass  transfer  and  spin  on  the  surface 
pressure  along  the  body. 
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Figure  24:  Case  2  Effects  of  mass  transfer  and  spin  on  the  heat 
transfer  along  the  body. 
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Figure  25:  Case  2  Effects  of  mass  transfer  and  spin  on  the  heat  transfer 
around  the  body  at  s  =  30.06. 
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Figure  26:  Case  2  Effects  of  mass  transfer  and/or  spin  on  the  heat 
transfer  around  the  body  at  s  *  30.06. 


11-233 


CM 


•  |  i  i  r 


Figure  28:  Case  2  Effects  of  mass  transfer  and  spin  on  the  streamwise 
sheat  around  the  body  at  s  -  30.06. 
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Figure  30:  Case  2  Effects  of  mass  transfer  and  spin  on  the  crossflow 
shear  around  the  body  at  s  *  30.06. 
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Figure  32:  Case  2  Effect  of  spin  on  the  crossflow  velocity  at 
s  =  6.40,  <f>  -  160  deg 
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Figure  33:  Case  2  Effect  of  mass  transfer  on  the  crossflow 
velocity  at  s  a  5.’85,  $  ■  160  deg. 
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